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The Square Kilometre Array (SKA) will transform our understanding of the role of the cold, 
atomic gas in galaxy evolution. The interstellar medium (ISM) is the repository of stellar ejecta 
and the birthsite of new stars and, hence, a key factor in the evolution of galaxies over cosmic 
time. Cold, diffuse, atomic clouds are a key component of the ISM, but so far this phase has been 
difficult to study, because its main tracer, the HI 21 cm line, does not constrain the basic physical 
information of the gas (e.g., temperature, density) well. 

The SKA opens up the opportunity to study this component of the ISM through a complementary 
tracer in the form of low-frequency (<350 MHz) carbon radio recombination lines (CRRL). These 
CRRLs provide a sensitive probe of the physical conditions in cold, diffuse clouds. The superb 
sensitivity, large field of view, frequency resolution and coverage of the SKA allows for efficient 
surveys of the sky, that will revolutionize the field of low-frequency recombination line studies. 
By observing these lines with the SKA we will be able determine the thermal balance, chemical 
enrichment, and ionization rate of the cold, atomic medium from degree-scales down to scales 
corresponding to individual clouds and filaments in our Galaxy, the Magellanic Clouds and be¬ 
yond. Furthermore, being sensitive only to the cold, atomic gas, observations of low-frequency 
CRRLs with the SKA will aid in disentangling the warm and cold constituents of the HI 21 cm 
emission. 
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1. The cold neutral medium 

The ISM is central to the evolution of galaxies as the birth site of new stars and the repos¬ 
itory of processed stellar ejecta. Star formation slowly consumes matter from the ISM and later 
returns it to the ISM via winds and explosions, enriched by the nucleosynthetic products of stellar 
interiors. Over cosmic time, this constant recycling and enrichment drives galaxy evolution. Stars 
also control the radiative energy budget of the ISM and its emission characteristics. Photons from 
massive stars with energies above 13.6 eV ionize hydrogen creating HII regions. Less energetic 
photons couple to the gas through photoelectrons from large molecules and small dust grains that 
heat atomic gas in PhotoDissociation regions (PDRs) surrounding HII regions and - on a much 
larger scale - in diffuse interstellar clouds. Through their winds and explosions, stars also stir up 
the ISM dynamically, sweeping up gas and forming large bubbles. This injection of mechanical 
energy into the ISM is a source of turbulent pressure, supporting the gas disk and the clouds therein 
against galactic- and self-gravity. Radiative and mechanical interaction of stars with their envi¬ 
ronment results in a number of distinct phases in the ISM; e.g., a phase of cold, neutral atomic 
hydrogen - diffuse clouds or Cold Neutral Medium (CNM) - and a phase of warm, neutral atomic 
hydrogen - the Warm Neutral Medium (WNM). 

Clearly, there is a complex feedback between stars and the ISM that determines the structure, 
composition, chemical evolution, and observational characteristics of the ISM in the Milky Way 
and other galaxies all the way back to the first galaxies. Understanding this interaction has to start 
with characterizing the phases of the ISM and linking these to stellar and galactic energy sources 
and physical processes that couple interstellar gas to the mechanical and radiative energy inputs 
from stars and the large scale motions in the galaxy. 

Over the last decades, our understanding of the ISM and its role in the evolution of galaxies has 
increased dramatically, but in leaps, driven, in particular, by the opening up of new observational 
windows by unique observatories. While the molecular medium has been well characterized and 
its role in star formation is well documented, analyzing and interpreting the atomic phase from the 
HI 21 cm line data in terms of the density and temperature has been challenging and has sofar only 
be done in a very limited number of cases (e.g. Heiles & Troland 2003). The main difficulty here 
is to accurately separate the warm and cold components that both contribute to the observed HI 
21 cm signal. A consequence of this is that the thermal balance of the cold, atomic ISM and its role 
in starformation is poorly understood. 

A promising, complementary method is to observe low-frequency CRRLs. These CRRLs have 
been found to be a good tracer of cold atomic gas (e.g. Payne et al. 1989, Anantharamaiah et al. 
1994, Kantharia & Anantharamaiah 2001, Anish-Roshi et al. 2001). Their observed properties, 
i.e. line strength and width, as a function of frequency provides a convenient method to determine 
the physical conditions of this gas. The unprecedented sensitivity and bandwidth of the SKA will 
revolutionize the field of low-frequency, i.e. < 350 MHz, CRRLs. In Sect. 2 we give a brief 
overview of the current status of the field. In Sect. 3 we discuss the models necessary to interprete 
CRRL observations. In Sect. 4 we describe how SKA-LOW will impact and transform the field 
of low-frequency CRRL sfudies. We focus our discussion here on two prime targets, the Galactic 
center and the Magellanic Clouds. In Sect.’s 5 and 6 we discuss the capabilities of SKA-LOW in 
phase 1 and 2 respectively, and finally we present a brief summary in Sect. 7. 
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2. Low-frequency CRRLs 

Radio recombination lines (RRL) result from atoms recombining with electrons in an ionised 
plasma and can broadly be separated in to two classes: discrete and diffuse (e.g. see Gordon & 
Sorochenko 2009 for a recent overview). Discrete RRLs from hydrogen (HRRL), helium (HeRRL) 
and carbon (CRRL) are typically observed above about 1 GHz. Here the RRL spectrum is domi¬ 
nated by HRRL and HeRRL lines that are found to be excellent tracers of the dense, warm, ionized 
gas (e.g. Palmer & Zuckerman 1966). Discrete RRLs are described in more detail by M. Thompson 
et al. (2014) and R. Beswick et al. (2014) in this volume. 

In this chapter we focus on the lower-frequency, i.e. <350 MHz, diffuse RRLs. At these 
frequencies we predominantly observe RRLs from carbon and hydrogen (e.g. Konovalenko & 
Sodin 1981, Anantharamaiah 1985.). HeRRLs have not been detected at these frequencies. Recent 
overviews of low-frequency RRLs and the instruments that are able to observe them are given 
in Peters et al. (2011) and Gordon & Sorochenko (2009). Below we will briefly summarize the 
current status of low-frequency RRL observations. 

2.1 Galactic CRRLs 

Large-scale studies: A summary of CRRL observations below 200 MHz is presented in Peters et 
al. (2011). Two large-scale Galactic surveys below 200 MHz were done at 76 MHz (Erickson et al. 
1995) and 34.5 MHz (Kantharia & Anantharamaiah 2001). These low-resolution (beam FWHM 
>2 deg) studies find that CRRL absorption is widespread throughout the inner Galactic plane with 
typical peak optical depths in the range 0.5 x 10^^ to 1 x 10^^ and typical line widths in the range 
10 to 50 km/s. HRRLs were not detected by these surveys. 

Two other low resolution (beam FWHM ~2 deg) surveys were conducted near 330 MHz 
by Anantharamaiah (1985) and Roshi et al. (2000). These surveys also find widespread RRF 
emission, from hydrogen and carbon, in the inner Galactic plane with peak optical depths in the 
range -0.5x10^^ to -2x10^^ for HRRF and -0.5x10^^ to -1x10^^ for CRRF (the minus sign 
indicates emission). The CRRFs are found to have typical linewidths in the range 15-30 km/s, 
about a factor 2 narrower than the corresponding HRRFs with typical linewidths 30-60 km/s. 

The broader HRRFs likely arise in the warm ionized medium (WIM) (e.g. Anantharamaiah 
1985, Heiles et al. 1996, Anish-Roshi et al. 2001, Kantharia et al. 2007). The CRRFs are more 
likely associated with the CNM as shown by their good correspondence with HI 21 cm absorption 
(e.g. Anish-Roshi et al. 2001, Anantharamaiah et al. 1994, Oonk et al. 2014). Theoretical models 
of CRRF emission are consistent with this picture (see Sect. 3). 

Kantharia & Anantharamaiah (2001) made an attempt at modelling the large-scale CRRF 
emission by combining the above mentioned surveys. They found that the observed CRRFs are 
consistent with arising in the CNM. However, their analysis was limited by the low resolution, 
poor spatial overlap and narrow frequency sampling of these surveys. Especially the unknown 
beam dilution due to the low resolution affected their investigation. Higher spatial resolution ob¬ 
servations are therefore essential to make further progress in this field. 

Small-scale studies: The only source that has been studied with sufficient spatial resolution, spec¬ 
tral coverage and sensitivity to allow for detailed CRRF modeling are the cold clouds along the line 
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of sight towards Cassiopeia A (Cas A). The CRRLs towards this souree ean reaeh a peak optieal 
depth of almost 10^^ (e.g. Kantharia et al. 1998). This is a faetor of a few brighter than seen on 
large seales in other regions, although this eould simply a eonsequenee of the low spatial resolution 
used in the studies of these other regions. The CRRL line width observed towards Cas A varies 
from about 4 km/s at 350 MHz to about 10 km/s at 50 MHz (Payne et al. 1989; Oonk et al. in 
prep.). 

Modeling shows that the CRRLs observed towards Cas A arise in the CNM (e.g. Payne et al. 
1989; Kantharia et al. 1998). Cas A is an example of a bright Galaetie souree that ean be used for 
detailed high-resolution CRRL ’pinhole’ studies of the CNM. Similarly extragalaetie sourees ean 
also be used to perform sueh pinhole studies, as shown by our reeent detention of CRRL absorption 
from the Milky Way against the bright extragalaetie souree Cygnus A (Oonk et al. 2014). 

With LOFAR we are eondueting a Galaetie survey to establish the physieal seales on whieh 
CRRLs are present. Some preliminary results are shown in Fig. 1. These and previous observations 
show that CRRLs are widespread in our Galaxy with typieal peak optieal depths of a few times 10~^ 
to a few times 10^^ and typieal line widths from a few to a few tens of km/s. 

2.2 Extragalaetie CRRL studies 

The first detention of extragalaetie low-frequeney CRRLs has reeently been aehieved with 
LOFAR. In Morabito et al. (2014) we report CRRL absorption from M82 near 60 MHz with a 
peak optieal depth of about 3x 10^^ and we find a good eorrespondenee with HI 21 em absorption. 
This first detention shows the potential that with inereasing sensitivity and bandwidth it is now 
becoming possible to detect CRRLs in other galaxies. In the upcoming years we will use LOFAR 
to carry out a flux limited, i.e. Fy >10 Jy at 160 MHz, survey for extragalaetie CRRLs. The 
primary aim of this survey will be to establish what type and fraction of bright extragalaetie radio 
sources show detectable CRRL signals and how the CRRLs compare to other properties of these 
systems. 

3. CRRL models & methods 

Low-frequency CRRLs originate in the partially ionized (Xg ~ 10~^), dense atomic phases 
of diffuse clouds and PDRs. This can be understood as singly ionized carbon has an ionization 
potential below that of hydrogen and this is therefore the dominant ionization state of carbon in the 
ISM. The observed properties of CRRLs, i.e. line strength and width, provide a convenient probe 
of the physical conditions of the gas. However, low-frequency CRRLs are not in local thermal 
equilibrium (LTE), and therefore non-LTE models are necessary to analyze the observations in 
terms of electron temperature, electron density and carbon column density (e.g. Shaver 1975; 
Brocklehurst & Salem 1977; Walmsley & Watson 1982; Ponomarev & Sorochenko 1992; Storey 
& Hummer 1995; Salgado et al. in prep.). A comprehensive literature exists on the topic of non- 
ETE models for low-frequency RRE emission and for further details we refer the interested reader 
to the above references and references therein. Here we show a few examples of such models in 
Eig. 2 and we briefly discuss how these models can be used to obtain the physical conditions of the 
CRRE emitting gas. 
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Electron temperature and density. The optical depth of CRRLs is affected by departure coef¬ 
ficients from LTE. For a given temperature Tg and density n^ these coefficients trace out a unique 
curve as function of frequency (or equivalently quantum number n), see Fig. 2. Therefore by mea¬ 
suring the CRRFs over a range of frequencies and comparing with models we can obtain Tg and n^ 
(e.g. Kantharia et al. 1998). Additional constraints on both of these parameters can be obtained by 
measuring how the pressure broadening affects the line width of the CRRFs (e.g. Kantharia et al. 
1998; Oonk et al. 2014). 

Thermal pressure: Once the electron temperature and density of the CRRF emitting gas has 
been obtained one can estimate the associated thermal pressure by invoking a carbon to hydrogen 
abundance [C/H]. The thermal pressure is then given by pth ~ngXTg/[C/H] (e.g. Payne et al. 1989, 
Kantharia et al. 1998). 

Carbon abundance: A systematic uncertainty in deriving hydrogen densities from CRRFs and 
thus thermal pressures is related to the carbon abundance [C/H]. Typically [C/H] has an uncertainty 
of about a factor 2 in our Galaxy. However, [C/H] can be determined by comparing, for the same 
line of sight, the carbon column density obtained from CRRFs with the cold hydrogen column 
density obtained from HI 21 cm absorption studies. 

WNM to CNM fraction: Fig. 2 shows that low-frequency CRRFs only trace cold, atomic 
gas, i.e. the CNM. This means that CRRFs can be used as a tracer of this phase in terms of 
temperature and distribution. To also obtain the hydrogen density it is necessary to invoke a carbon 
to hydrogen abundance [C/H]. A Galactic [C/H] grid can be created by comparing CRRF and HI 
21 cm absorption measurements (e.g. this work; McClure-Griffiths et al. 2014 in this volume). If 
[C/H] does not vary strongly between the HI 21 cm absorption measurements then the CRRFs can 
be used to also determine the CNM hydrogen density between absorption measurements. 

Hydrogen ionization: Hydrogen in diffuse clouds can only be ionized through cosmic rays. 
Therefore the HRRF to CRRF ratio, if they measure the same gas, is a measure of the cosmic ray 
ionization rate (e.g. Sorochenko & Smirnov 1987). Recently, using observations from the WSRT 
at 350 MHz and FOFAR at 220 MHz, we have found a weak, narrow HRRF component in the 
cold clouds towards Cas A. These HRRFs have the same width and velocity as the corresponding 
CRRFs indicating that the arise from the same gas and thus can be used to determine the hydrogen 
ionization rate, see Fig 3. 

4. Square Kilometre Array 

The properties of the low-frequency CRRFs discussed in Sect. 2 & 3 show that these originate 
in the CNM where they are associated with cold, diffuse clouds and PDRs. Here we propose to 
use SKA-FOW to perform a survey of Carbon Radio Recombination Fines (CRRFs) originating 
from cold atomic gas in diffuse interstellar clouds and in PDRs associated with regions of massive 
star formation in the Milky Way as well as with AGN and starburst activity in other galaxies. This 
CRRF survey can help answer a number of important open questions related to the ISM, such as; 

(1) What is the morphology of CRRL emitting regions and how do they relate with the neutral, 
molecular, star forming and hot gas, as well as the large-scale structure of the Milky Way? 
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(2) What is the thermal, pressure balance for cold, diffuse clouds? 

(3) What is the C abundance in the diffuse gas and how does this vary with galactic radius? 

(4) What is the ionization rate for the cold ISM? 

(5) What are the physical conditions of the cold, atomic gas in other galaxies and how do these 
compare to other properties like jets, winds, and stellar, molecular & dust mass? 

Instead of attempting to summarize here all of the possible science topics that can be pursued 
with CRRLs and the many interesting regions that can be investigated we have chosen to limit our 
discussion in this section to two prime targets in the southern hemisphere, the Galactic center and 
the Magellanic clouds. In the Sect.’s 5 and 6 we will then discuss the technical details for studying 
these objects with SKA-LOW in phase 1 and phase 2. 

SKA & Galactic Center. The center of the Milky Way shares many properties with other, often 
more spectacular nuclei of galaxies. Indeed, at a distance of only 8 kpc, the nucleus of our galaxy 
is a unique laboratory for studies of global phenomena - e.g., accretion, mass flow, star formation - 
and the physical and chemical conditions and the micro-physics processes at work in the interstellar 
medium in the nuclei of galaxies. The few parsecs surrounding the massive black hole in the center 
of our Milky Way contain a dense and luminous star cluster, as well as various components of 
neutral, ionized and hot gas. 

The central parsec region contains mostly ionized gas (the HII region Sgr A West) and hot, X- 
ray emitting gas. This low density, ionized ’central cavity’ is pervaded by a set of orbiting ionized 
gas filaments, which in turn are surrounded by orbiting, dense molecular cloud streamers at ~ 1.5 to 
4 pc (the circum-nuclear disk). The outer edge of the circum nuclear disk is bordered by the young 
supernova remnant, Sgr A East, and surrounded by a number of dense molecular clouds on a scale 
of 5 to 100 pc (see the reviews by Genzel et al. 2010; Morris & Serabyn 1996, and references 
therein). This central molecular zone is embedded in a more extensive HI nuclear disk. The high 
surface density and total mass content likely reflect inflow from larger radii. 

Because of the high pressures and strong tidal shearing effects, the average properties of the 
molecular clouds in the Galactic Center differ significantly from those in the outer Galaxy with 
elevated gas temperatures, high average densities and strong turbulent velocity fields. Magnetic 
fields with mG field strength have been reported, that - at least locally - control the dynamics 
of the gas. A number of gas heating mechanisms are important in this nuclear environment in¬ 
cluding far-UV radiation emitted by the OB stars near the center responsible for the HII regions 
and the surrounding photo-dissociation regions at the edges of the cavities blown by the hot stars. 
Clouds compressed at the edges of the expanding bubbles are shock heated. The dense gas phase 
in shielded cloud cores is heated by turbulent dissipation and magnetic viscous heating. The inter¬ 
relationship of these structures and the processes that drive them will be an important testbed for 
our general understanding of what drives the evolution of galactic nuclei. 

Low-frequency CRRLs can provide a unique information on the central region of our galaxy 
specifically on the physical conditions in and the dynamics of the dense photodissociation regions 


6 




Low-frequency RRLs with the SKA 


J. B. R. Oonk 


that separate the ionized and moleeular gas (Wyrowski et al. 1997, 2000; Natta et al. 1994). On a 
larger scale, the carbon radio recombination lines can probe the interrelationship of the molecular 
gas and the wider atomic H layer. 


The Magellanic Clouds. The Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) 
are irregular galaxies in orbits around the Milky Way that have experienced several tidal encounters 
with each other and with the Galaxy. The Large Magellanic Cloud contains a bright stellar bar 
which together with regions of star formation and the patchy dust obscuration dominates its optical 
appearance. In HI, on scales of 15 to 500 pc, the morphology is dominated by HI filaments with 
numerous shells and holes. The turbulent and fractal nature of the ISM on these scales is the result 
of dynamical feedback into the ISM from star formation processes. On larger scales, (200 pc to 
10 kpc), the HI distribution is very patchy. The reduced metallicity implies lower dust shielding 
and higher photodissociation rates. Hence, the Magellanic Clouds are a more hostile environment 
for molecular gas. Indeed, while - because of self-shielding - much of the hydrogen may be in 
molecular form in the clouds, carbon may be largely in the form of C+. In general, this so-called 
CO-dark gas may be an important reservoir of gas in low metallicity environments. 

The LMC is home to the most luminous and best known massive star-forming region in the 
Local Group, the Tarantula nebula, powered by the super star cluster 30 Dor, which contains some 
1000 bright O stars. The Magellanic Clouds serve as excellent laboratories for the study of the 
process of star formation in an environment that, in many ways (such as metallicity, dust content, 
and star formation rate), resembles the extreme conditions of the early universe and distant star¬ 
forming regions. 

The LMC and SMC are also excellent astrophysical laboratories for studying the life cycle of 
the interstellar medium (ISM). Located at 50 and 60 kpc, their proximity permits detailed studies 
of resolved ISM clouds and their relation to stellar populations on global scales, which in the Milky 
Way are obfuscated by line-of-sight crowding and distance ambiguity. The SMC and LMC (and 
Milky Way) span an interesting range in their global properties (metallicity, molecular gas content, 
and gas-to-dust ratio) that make them interesting testbeds for galaxy evolution studies. In terms of 
galaxy evolution, the LMC-SMC pair is well suited for studying how the agents of evolution, the 
ISM and stars, operate as a whole in two galaxies that are tidally interacting with each other and 
with the Milky Way. 

Low frequency carbon recombination lines will provide a unique tool specifically geared to¬ 
wards tracing CO-dark gas and its relationship to, on one hand, HI clouds and molecular gas and, 
on the other hand, the regions of active star formation. In addition, CRRLs can be used to study the 
physical conditions in the photodissociation regions associated with regions of massive star forma¬ 
tion. In particular, the 30 Dor region will allow a study of the characteristics of PDRs in a super star 
cluster environment. Finally, we like to highlight that the many bright Hll regions and supernova 
remnants will provide convenient bright background sources for detailed studies of the carbon and 
hydrogen recombinations that can probe the cosmic ray ionization rate on a galaxy-wide scale and 
relate them to their likely sources. 
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5. SKA Phase 1 

The technical details of SKA-LOW in phase 1, SKAl-LOW, are described in the SKAl Imag¬ 
ing Science Performance document (hereafter SKAl-ISP) by R. Braun. SKAl-LOW provides an 
increase in sensitivity, relative to its precursors and pathfinders, by up to 2 orders of magnitude. 
The proposed configuration offers a dense core, with a filling factor close to unity. Outside the core 
the stations are arranged in three outwards extending spiral arms. The core dominates the angular 
resolution versus sensitivity and the sweetspot is in range 1-5 arcmin at 140 MHz. SKAl-LOW 
will have a large bandwidth from 50-350 MHz where up to 243 a-transition CRRLs (n=266-508) 
can be observed simultaneously. Comparing with the models in Fig. 2 shows that this range is well 
suited to study the low density, low temperature CRRL emitting gas. 

Galactic survey. The overall goal of a Galactic CRRL survey with SKAl-LOW is to determine 
the physical conditions of the cold, diffuse clouds and their role in the process of starformation in 
the inner Galactic plane (1<50 deg, b<l5l deg). In Sect. 2 we showed that, at frequencies below 
350 MHz, on degree-scales Galactic CRRLs have typical line widths in the range from 10 to 50 of 
km/s and typical peak optical depths in the range of a few times 10^^ to a few times 10^^. 

Table H] shows that upon stacking the CRRLs in sets of 9 lines that the SKAl-LOW sensitivity 
in 8 hours with 1 kHz resolution is sufficient to detect CRRLs down to a peak optical depth level 
of ~10^^ throughout the entire inner Galactic plane on scales from 1.5 arcmin at 350 MHz to 10 
arcmin at 50 MHz. The arcmin-scale resolution is well suited to study the cold, diffuse clouds 
in the Milky Way and the peak optical depth limit corresponds to a HI column density of about 
IxlO^o cm-2. 

Stacking CRRLs in sets of 9 does not affect our science goals as the lines are spaced very 
closely. In this way we obtain 27 stacked CRRL a transition measurements for each line of sight 
across the full frequency range. In Fig. 2 we give an example of what such a measurement with 
SKAl-LOW simulated noise could look like. The poorest sensitivity, in units of optical depth, is 
obtained at 350 MHz. However, by convolving the higher frequency maps to the spatial resolution 
of the lowest frequency maps at 50 MHz the noise will be similar at all frequencies. By comparing 
the SKAl-LOW noise levels with our models this means that we can determine the electron tem¬ 
perature and density on arcmin-scales in diffuse clouds to better than a factor 2 down to a limiting 
HI column density of about 1 x 10^° cm^^. 

Additional constraints on the electron temperature and density of the cold gas can be obtained 
from the CRRL line width (e.g. Kantharia et al. 1998). However, the currently planned 1 kHz 
channel width in phase 1 is not sufficient to resolve the narrowest CRRLs measured towards Cas A 
(e.g. Fig. 3; Payne et al. 1989) and Cygnus A (Oonk et al. 2014). To resolve these lines would 
require a channel width of about 0.3 kHz. 

A 500 deg^ survey of the inner Galactic plane requires 13 pointings at 50 MHz and 500 point¬ 
ings at 350 MHz. For 8 hr pointings this would take about 10^ hrs at 50 MHz and 4x10^ hrs at 
350 MHz and is therefore limited by the small field of view at 350 MHz. If a beamformer becomes 
part of the SKAl-LOW design this time requirement can be reduced by scanning the sky in smaller 
frequency chunks with more beams. 

One of the prime targets of a Galactic CRRL survey with SKAl-LOW is the Galactic center. 
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At the distance of the Galactic center (8 kpc) the 1.5 arcmin resolution at 350 MHz means that we 
can resolve CRRL emitting structures down to 4 pc at 350 MHz. To study the physical properties 
of these structures we are limited to the 10 arcmin resolution at 50 MHz corresponding to about 
25 pc. This means that in phase 1 we can study the interaction of the molecular gas with the wider 
atomic H layers. Smaller scale structures such as the PDRs and the HI nuclear disk require the 
capabilities of SKA-LOW in phase 2. 

To derive the carbon abundance and convert our carbon column densities to hydrogen column 
densities we aim to combine our CRRL measurements with the planned Galactic HI absorption grid 
survey, see McClure-Griffiths et al. (2014) in this volume. We estimate that this can be done on 
average in the Galactic plane for about 5 extragalactic background sources per deg^ (e.g. Wilman 
et al. 2008; Riley et al. 1999). This number should be considered a lower limit for measuring 
the carbon abundance in the galactic center given the many bright Galactic radio sources that are 
located in this region and which can also be used for this purpose. 

The models in Fig. 2 and the Gas A spectra in Fig. 3 show that with with the spatial resolution 
and sensitivity of SKA 1-LOW it should be possible to detect the narrow HRRL component from 
diffuse clouds and separate it from the broader HRRL components associated with the WIM. The 
models show that this is best done at higher frequency range of SKAl-LOW. In combination with 
the corresponding CRRLs this will allow us to determine the hydrogen ionization rate due to cos¬ 
mic rays in these clouds. Presently, this rate in the CNM has only been measured through molecular 
observations in the infrared or submillimeter (e.g. Indriolo et al. 2007; Shaw et al. 2008). How¬ 
ever, these observations have been limited to a small number of bright sources and the analysis 
is hampered by the complex molecular physics involved. Cosmic ray ionization, especially in the 
Galactic center, may play an important role in heating of diffuse clouds. 

Extragalactic survey. The overall goal of an extragalactic CRRL survey with SKAl-LOW is 
to determine the presence of CRRLs in other galaxies and use these lines to trace the physical 
conditions of the CNM in these objects. Our CRRL models show that we can expect a peak optical 
depth of about 10~^ for a HI column density of about 10^^ cm^^ (Salgado et al. in prep.). This 
is consistent with our recent detection of CRRLs in M82 (Morabito et al. 2014). Furthermore 
the known frequency spacing of RRLs can be used to determine redshifts for extragalactic objects. 
This is particularly useful for those optically faint objects where it is difficult to identify the HI 
21 cm line. 

Table || shows that upon stacking the CRRLs in sets of 9 lines that the limiting point source 
sensitivity is about 0.4 Jy at 160 MHz for an 8 hr observation with SKAl-LOW and a resolution of 
about 3 arcmin. A 2% sr survey to this depth would enable a CRRL search in up to 10^ galaxies, 
spanning a redshift range from z=0 to 5 (Wilman et al. 2008). This sample size would represent an 
increase of almost 3 orders magnitude relative to our planned LOFAR survey. 

The time requirement for such a 2% sr survey varies between 4x10^ hrs at 50 MHz and a 
staggering 1.6x10^ hrs at 350 MHz. At the high frequencies this is clearly too long to be real¬ 
istic. If a beamformer becomes part of the SKAl-LOW design then these time requirements can 
be reduced significantly. Alternatively, one could first perform a blind CRRL detection survey 
near 100 MHz and then only follow-up the detected galaxies with wider frequency coverage using 
pointed observations. 
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Of particular interest in an extragalactic SKAl-LOW survey are nearby galaxies. We estimate 
that for 8 hr integrations about 400 nearby galaxies are bright enough to be detected CRRLs and 
that in about 100 of these the CRRLs can be spatially resolved (e.g. Israel 1990, Negrello et al. 
2013). A pointed survey for nearby galaxies over the full frequency range would therefore take 
about 4x 10^ hrs. 

Two prime targets in a nearby galaxy survey are the Magellanic Clouds. They contain large 
quantities of cold, atomic gas (e.g. Stanimirovic et al. 1999; Dickey et al. 2000) and have a 
408 MHz background radio continuum temperature of about 30-120 K with an approximate ^ 
scaling to lower frequencies (Haslam et al. 1982; Israel et al. 2010). Table then shows that in 
an 8 hr single pointing we can reach peak optical depths of few times 10^^ at 350 MHz to a few 
times 10^“^ at 60 MHz. This is sufficient to detect CRRLs for hydrogen column densities down a 
few times 10^° cm^^ on 10 arcmin scales. 

At the distance of the Magellanic Clouds a 10 arcmin scale corresponds to 200 pc. This means 
that we will be able to resolve the CRRLs associated with larger scale patchy HI distribution and 
HI supershells. The smaller scale HI shells and PDRs require the capabilities of SKA-LOW in 
phase 2. Just as for the inner Galactic plane survey we estimate there are on average 5 extragalactic 
background sources per deg^ for which HI 21 cm absorption and CRRL measurements can be 
compared to derive the carbon abundance in these systems. 

5.1 SKA Phase 1: Early science 

For early science (ES) it is expected that SKAl-LOW will have about half of its final sensifiv- 
ify. The CRRL science fhaf can be done during ES depends fhe array configuralion. If fhe filling 
factor of core of fhe array during ES remains similar fo fhaf of fhe final array configuralion Ihen 
fhe inner Galaclic plane survey, as ouflined above, can be carried ouf. Given fhe high HI column 
densifies and high background confinuum radiation fhe survey can fhen sfarf al fhe Galaclic cenler 
during ES and progress oulwards in Galactic longilude as fhe array builds up to ils full phase 1 
sensilivily. 

Eor an exlragalacfic CRRE survey fhe lower sensilivily during ES has differenl effecls for 
differenl source populations. Eor fhe brighf, dislanf ER I and ER II radio galaxies fhe number of 
objecls fhaf can be searched for CRREs down lo a peak oplical deplh of 10^^ decreases wilh aboul 
a faclor of Iwo (Wilman el al. 2008). Eor nearby galaxies Ihe effecl is more severe. We estimate a 
decrease by aboul a factor 4, such lhal al mosl 100 nearby galaxies are brighf enough to be searched 
for CRRLs and lhal only aboul 20 of Ihese are sufficienlly brighf and large lhal Ihe CRRLs wilhin 
can be spatially resolved (Negrello el al. 2013, Israel el al. 1990). 

Eor Ihe Magellanic Clouds Ihe lower sensilivily during early science will primarily affecl Ihe 
limiting column density. Therefore for a fixed observing time of 8 hrs per pointing we will slarl to 
miss Ihe more diffuse clouds slarling al a HI column density of aboul 5 x lO^*’ cm^^. 

6. SKA Phase 2 

In phase 2 SKA-EOW is expected to have 4 times Ihe sensitivity of SKAl-EOW over a wide 
range of angular scales and up to 20 times Ihe angular resolution. Eor bolh a Galactic and an 
exlragalacfic survey Ihis increase in resolution is very imporlanl. In Ihe Galactic center we will be 
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able to spatially resolve the HI nuelear disk and eonneet its kinematieal and thermal properties to 
that of wider inflow of eold gas. Similar this inerease in resolution is essential to spatially resolve 
PDRs and HI shells, as traeers of stellar feedbaek, over a wide range of environments in both our 
own Galaxy and the Magellanie Clouds. 

For a larger extragalaetie survey the inereased sensitivity and resolution of SKA-LOW in phase 
2 would allow us spatially resolve CRRLs in up to 3x 10^ nearby galaxies and seareh for CRRLs 
in up to half a million >0.1 Jy radio sourees spanning a range in redshift from z=0 to 5. This 
will provide important input on how the physieal properties of the CNM evolve with time and 
environment. If redshift information is available then the average CNM properties for well-defined 
samples of even fainter radio sourees ean be studied by employing staeking teehniques. 

7. Summary 

SKA observations of low-frequeney CRRLs, in eombination with HI 21 em absorption mea¬ 
surements, will provide a eomprehensive and unpreeedented inventory of the low-density, eold 
ISM in our Galaxy, the Magellanie Clouds and beyond. This will enable us to quantify its role in 
the overall pressure, mass and energy balanee of the interstellar medium and how it relates to the 
cyeles of star formation and death. 

In addition, for studies of the Epoeh of Reionization and the HI dark ages it is important to 
quantify the distribution of the Milky Way CNM, in terms of strength and veloeity, as a potential 
foreground signal. Deep, high speetral resolution, low-frequeney RRL observations with SKA- 
LOW ean help provide this quantifieation and by doing this with the same instrument at the same 
frequeneies one ean minimize the uneertainties due to instrumental systematies and extrapolation 
from higher to lower frequeneies. 

In order to speetrally resolve CRRLs assoeiated with the eoldest, atomie gas struetures on 
small seales a ehannel width of about 0.3 kHz will be required, instead of the 1 kHz that is eurrently 
planned in the SKAl design. Another important element that eould signifieantly speed up large 
low-frequeney surveys is a beam former. The ehange in field of view with frequeney for SKA- 
LOW means these surveys ean be done mueh more effieiently by trading bandwidth against beams. 
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Frequency 

^ T, rms , ext end ed 

peak,extended 

,rms,point 

F V,point 

[MHz] 

[mK] 

x(10-4) 

[mJy] 

[Jy] 

60 

980 

(0.5-0.05) 

1.9 

1.9 

160 

140 

(1.0-0.10) 

0.4 

0.4 

220 

140 

(2.1-0.21) 

0.2 

0.2 

350 

154 

(7.5-0.75) 

0.2 

0.2 


Table 1: CRRL sensitivity with SKAl-LOW. Column (1) Observing frequency. Column (2) Core SKAl- 
LOW 5(7 spectral sensitivity for extended sources, in units of mK, for an effective baseline Bg//=2 km with 
1 kHz channels, 8 hr integration and stacking 9 lines. Column (3) Corresponding 5(7 peak optical depth range 
for extended sources. These ranges are calculated for the inner Galactic plane with continuum temperatures 
between 140 and 1400 K at 408 MHz (Haslam et al. 1982) and using a scaling (e.g. Rogers et al. 2008). 
Column (4) Full SKAl-LOW 5(7 point source spectral sensitivity, in units of flux, for 30 km/s channel 
width, 8 hr integration, 200 arcsec FWHM beam and stacking 9 lines. Column (5) Corresponding limiting 
point source continuum flux Fy,/™ to obtain a 5(7 CRRL detection with Comparing to the S3 

simulation (Wilman et al. 2008) and the 7C survey (Riley et al. 1999) shows that on average there are about 
5 extragalactic sources per deg^ brighter than these limiting fluxes. For both the extended and point source 
cases SKAl-LOW is typically about a factor 10 more sensitive than LOFAR at 160 MHz and a factor 100 
more sensitive than LOFAR at 60 MHz. The above sensitivities are derived from the SKA-ISP document. 


Cyg X: Integrated CRRL optical depth [km/s] 
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Cos A: integrated CRRL optical depth [km/s] 
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Figure 1: Observations of CRRL absorption in the Milky Way near 50 MHz (Oonk et al. in prep.). (Left) 
LOFAR 4 hr observation of degree-scale CRRLs in the Cygnus X region. Shown is a map of the integrated 
CRRL absorption in units of optical depth with radio continuum contours overlaid in white. The map shows 
structure down to the 1 degree limit of this observation. (Right) LOFAR 8 hr observation of arcmin-scale 
CRRLs along the line of sight to Cas A. Shown is a map of the integrated CRRL absorption in units of 
optical depth with radio continuum contours overlaid in white. The CRRL absorption is resolved in several 
filamentary structures that agree well with the corresponding HI 21 cm absorption (e.g. Bieging et al. 1991). 
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Figure 2: RRLs and the diffuse neutral ISM. Model results are given for the CNM (n//=50 cm 
ng=0.05 cm-3, T<,= 100 K) and WNM (n/r=0.1 ne=0.01 cm-^ 1^=10000 K). For both the CNM and 

the WNM a hydrogen column density N(HI)=10^* cm^^ is used.fLe/fj Integrated optical depth vs. quantum 
number for CRRL (solid) and HRRL (dashed) for the CNM (black) and WNM (green). The WNM optical 
depth is negligible and as such the HRRL and CRRL curves for the WNM overlap on this plot. (Right) 
The CRRL curve for the CNM model with error bars simulated for a 8 hr single pointing core SKA 1-LOW 
observation in the inner Galactic plane with a continuum temperature of 140 K at 408 MHz. The integrated 
optical depth depends linearly on the column density. In red and blue we overplot two lower temperature 
models that are scaled to reproduce the optical depth for the black reference model at n=390. 


Cassiopeia A (1.0 km/s chonnels) 



Cas A (HBA-HIGH) 



Figure 3: LOFAR and WSRT observations of CRRL and HRRL emission in the cold clouds towards Cas A. 
(Left) WSRT 12 hr 350 MHz observation. The RRL spectrum showing matching narrow HRRL and CRRL 
emission. For both the WSRT and LOFAR spectrum the velocity scale has been centered on the local 
standard rest for CRRLs. on this scale the corresponding HRRLs are offset by h- 150 km/s. (Right) LOFAR 
3 hr HBA-HIGH observation at 220 MHz. This spectrum was obtained using only the LOFAR core stations. 
The HRRLs in both the LOFAR and WSRT spectra are likely associated with hydrogen that has been ionized 
by cosmic rays and the HRRL/CRRL ratio implies a cosmic ray hydrogen ionization rate ~10^*® s^* 
(Oonk et al. in prep.). 
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